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1 Introduction 

There have recently been numerous conjectures regarding spectacular 
multi-TeV electroweak phenomena, characterized by the production of many 
electroweak gauge and Higgs bosons, which may be accessible to the next 
generation of colliders. Though such phenomena were initially studied in 
the context of topologically induced baryon plus lepton number violating 
interactions in the Standard Model 0, it was subsequently realized that 
they are a general feature of large order tree-level processes, where dramatic 
enhancements in production cross sections may be attributed to phase space 
and combinatorial factors 0. 

Yet another proposal for multiple weak gauge boson production was put 
forward recently by Ringwald et al.^ |J under the name of geometrical 
flavour interactions. They propose an analogy between multiparticle pro- 
duction in Quantum Chromodynamics (QCD), which occurs above energy 
scales of a few GeV, and possible multiparticle production in weak inter- 
actions above the TeV scale, where weak gauge boson masses may be ne- 
glected. More precisely, they relate QCD to a prototypical confined phase 
of weak SU(2) and, using scaling arguments, extrapolate the analogy to the 
Higgs phase. As a consequence they allege that above parton-parton center of 
mass energies of 2 to 20 TeV, the non-perturbative production of 0(a~ 1 ) ~ 30 
electroweak bosons could occur with cross sections of 0(1 nb - 10 //b). Ir- 
respective of the appropriateness of their scaling arguments, Ringwald et al. 
qualitatively suggest that, as far as weak interactions are concerned, collid- 
ing quarks and leptons may be thought of as being surrounded by clouds of 
virtual weak bosons which eventually result in "black disc" cross sections, 
and hence the name of geometrical flavour interactions. 

We remark that Ringwald et a/.'s qualitative view of high energy weak 
interactions is reminiscent of the QCD parton model in which an energetic 
proton may be viewed as a collection of quarks and gluons which form a 
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cloud of transverse size 0(Aqq D ). In contrast to the approach of Refs. || |], 
we investigate in this paper the possibility of high energy weak multiparticle 
production in the context of perturbation theory using techniques similar 
to those used for multiparticle production in QCD. We take two successively 
refined viewpoints. First, we adopt a weak parton model in which weak gauge 
bosons assume the role played by gluons in QCD and look for the failure of the 
lowest order approximations to the corresponding weak parton distribution 
functions which we interpret as signalling the increased importance of states 
with many gauge bosons. Second, we perform exact calculations of selected 
2 — > 2, 2 — ► 3, and 2 — > 4 weak processes to investigate particle production 
at mass scales inaccessible to the parton model. For illustrative purposes 
we concentrate on contributions to the v^P^ total cross section, a^J, chosen 
because it is naively free from QCD and QED complications. 

Our goal is to quantitatively explore the energy range in which multipar- 
ticle final states become a dominant feature of weak interactions. Since we 
limit ourselves to the first few orders of perturbation theory, we will only 
consider final states with up to three or four particles (as opposed to the 
production of 0(a~ x ) particles discussed in Refs. [JJ, M). Furthermore, we 
are interested in the perturbative decomposition of cr t ^ into its various mul- 
tiparticle channels. The more global question of the perturbative growth of 
total (inclusive) weak cross sections with \/i has been considered recently 
by Pindor, Raczka and Wetterich in Ref. J5] in which they solve an integral 
equation suggested by Fadin, Kuraev and Lipatov|| [7|. They find slowly 
growing weak total cross sections (i.e., 0(100 pb) for y/s < 10 15 GeV) but, 
due to the nature of their methods, they cannot decompose their results into 
contributions with definite numbers of particles in the final state. 

Using a variety of techniques we will argue that weak multiparticle pro- 
duction becomes a dominant characteristic of weak interactions at energy 
scales in the range 10 7 — 10 9 GeV which is firmly beyond the reach of any 
planned accelerator. Kane and Scanio M have pointed out that cross sections 



for exclusive processes such as e + e _ — > fi + fi~ can be surpassed by higher or- 
der processes (such as e + e~ — ► ^ e z/ e /i + /z~ through W + W~ fusion) at center 
of mass energies of O(10 3 GeV). However this dominance is only relative 
to a rapidly falling s— channel Born contribution. In contrast, lowest order 
VePfx scattering proceeds through t— channel exchanges which are analogous 
to dominant t— channel processes in QCD; at higher orders multiparticle 
production arises from very forward (i.e., small angle) inelastic reactions. 
We shall not consider the possibility of a strongly interacting Higgs sec- 
tor which could, by itself, dominate the characteristics of high energy weak 
interactions Q. 

The paper is organized as follows. In section 2 we review the QCD parton 
model approach for describing perturbative contributions to the antiproton- 
proton total cross section a^ tgl . In section 3 we consider the lowest order 
distribution functions of the weak parton model. We use them to determine 
contributions to cr^j through gauge boson fusion and investigate the energy 
scale at they become unreliable. In section 4 we present results of calculations 
beyond the weak parton model and argue that the dominance of multiparticle 
states at energy scales of 10 7 — 10 9 GeV is more than an artifact of large 
logarithms in finite order calculations. 

2 QCD Parton Model and cT t p p tal 

Let us abstract those aspects of the perturbative QCD parton model of 
relevance to high energy weak interactions. Underlying this ambition is the 
operating premise, which is not without controversy, that multiparticle pro- 
duction and the growth of crfotai with center of mass energy may be deduced 
from the QCD parton model. We briefly review the support and criticisms 
of the QCD parton model view of cr{^ al to set the stage for the weak parton 
model in the next section. 



Though an old idea|IIJ, perturbative contributions to rising hadronic 
cross sections gained prominence with the UAl observation of minijets[|TT] 



and have subsequently led to many theoretical studies of the perturbative 
component of cr^ al [|12| [H|. However, realistically the detailed successes of 
a perturbative parton model at momentum transfers at O(Aqcd) must be 
viewed with a degree of skepticism since there is little a priori justification 
for a perturbative treatment in that regime. Indeed, descriptions of c^ tal in 
terms of the QCD pomeron have also been shown to accommodate existing 
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As an approximation to cr^ al in the perturbative QCD parton model one 

d<ji 



^t P o P tai(s) - J2 L dp 2 T dxi dx 2 fi/ P {xi, Q 2 ) fj/p(x 2 , Q 2 ) -rf(Q 2 ) 

~f Jp 2 T >l GeV 2 J dp? 

(1) 

where Q 2 = X1X2S, fi/ p {x) is the number density of parton species i carrying 
a proton momentum fraction x (i.e., a parton distribution function) and 
aij(Q 2 ) is the subprocess cross section for colliding partons of species i and 
j. (It is actually more appropriate to evaluate the distribution functions at a 
scale 0(p1) but such details will not concern us.) The sum extends over all 
possible quark and gluon species. The ad hoc infrared cutoff p\ > 1 GeV 2 
avoids the non-perturbative growth of a s {p1 — > Aq CD ). 

In the parton framework the perturbative behaviour of o~f^ tal (s) is accom- 
modated as follows. As the center of mass energy ^ increases, successively 
smaller x values become kinematically accessible. For example, the small-x 
behavior of the gluon distribution function f g / p (x, Q 2 ) is thought to behave as 
~ x 5 (where 5 < 0) so that gluon-gluon collisions give large contributions to 



the total cross section|15|. The cooperation of these kinematic and dynamic 
effects may be illuminated by a toy calculation in which one considers only 
gluon-gluon collisions in Eq. [j], and takes f g / p (x) ~ 1/x and a gg ~ constant, 
which gives a tot ~ In s. The enhanced small-x behavior of parton densities 



inside hadrons suggests the picture of a high energy proton as a cloud of 
small-x partons. Indeed, using more realistic distribution functions and hard 
cross sections one can reproduce the Tevatron result for af^ al — 70 — 80 mb 
at -y/J = 1.8 TeV |16|| . The ability to mimic this result using perturbative 
techniques is based upon the assumption that the total cross section is dom- 
inated by semi-hard collisions; the infrared cut-off is then interpreted as the 
scale at which semi-hard processes become relevant. In the following section 
we exploit this apparent success of perturbative QCD and apply similar tech- 
niques to describe high energy weak interactions where weak gauge bosons 
assume the role played by gluons in QCD. 

3 The Weak Parton Model 

3.1 Factorization and Distribution Functions 

In QED there are already examples of rising cross sections which may 
be understood in the context of a parton picture. Among them is the two- 
photon production of a fermion pair (//) through e + e~ — ► e + e~ ff which 
has a lowest order cross section f]17 
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.(e+e-^eW/H^ln^jln^j, (2) 

which may be obtained by convoluting the probability of finding quasi-real 
photons inside electrons with the subprocess cross section for fermion pair 
production from two on-shell photons. This is the Weiszacker- Williams ef- 
fective photon approximation JI§|| , which is in good agreement with exact 
calculations |19|]. The parton model approach to QED has also been success- 



fully applied to describe other higher order electromagnetic processes such 
as corrections to the Z boson line shape |2^ . 



The enhanced small- a; behaviour of the gluon distribution function in 
QCD and the large QED logarithms in Eq. have counterparts which we 



wish to exploit by considering a parton model for the weak interactions. By 
analogy with QCD, we consider writing contributions to the cr^j in the 
factorized form 

°totai( s ) = Yl J dx ± dx * fi/»*( x u Q 2 ) fj/vpfa* Q 2 ) °ij(Q 2 ) ( 3 ) 

where Q 2 = X1X2S. In principle, the sum in Eq. [| extends over all weak 
parton species found "inside" neutrinos (e.g., W, Z, v, ■ ■ •). The weak par- 
ton distribution functions and the on-shell parton-parton subprocess cross 
sections have interpretations similar to those in QCD. At the level of this 
discussion we emphasize that the factorization of u^tai * n Eq- is performed 
strictly by analogy with QCD and as such is an assumption without rigorous 
justification^]. As in QCD, the scale at which the distribution functions are 



evaluated is not actually X1X2S but rather is determined by the subprocess 
under consideration. 

In contrast to the QCD factorization, the integral in Eq. [^ does not have 
an arbitrary infrared cut-off imposed by the dynamics of the theory. Instead, 
there are implicit kinematic constraints of the form 

(m. + m,)' 

s 

as required by the assumption that the weak partons participating in the 
hard interaction are on-shell with masses m$ and rrtj. In principle, similar 
constraints are present in QCD but are usually irrelevant since they are 
superseded by the more restrictive constraint on p 2 -, associated the onset of 
non-perturbative phenomena. 

The constraints of Eq. f| concern us for the following reasons. Drawing 
from our experience with QCD and QED, we anticipate that the small-x 
enhancement of the distribution functions of weak gauge bosons is a promis- 
ing place to look for cross section enhancements. However due to Eq. f| the 
weak parton model can only address gauge boson subprocesses above energy 
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scales of O(Mjy). It is uncertain beforehand whether the bulk of the mul- 
tiparticle production will occur above or below this scale. While multiple 
gauge boson production certainly requires energies of O(Mw), the produc- 
tion of single gauge bosons or light quarks and leptons can occur below this 
scale. For weak interactions we can, unlike QCD, always perform an exact 
perturbative calculation (and we do) to check the sensibility of our assump- 
tions. However we initially follow the more intuitive route of the weak parton 
model to elucidate the operative mechanisms. 

The weak parton distribution functions in Eq. |3| obey coupled evolution 



equations analogous to the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi|23[] equa- 



tions in QCD and the Lipatov equations in QED|24|]. Due to the rela- 



tively small coupling constants involved, the evolution equations for QED 
and weak interactions can be solved iteratively in terms of parton splitting 
functions (see Appendix A) starting from well defined input distributions 
(e.g., fv/ u (x,Q 2 = 0) = 5(1 — x) ). In contrast, the strong QCD coupling 
constant forbids such an expansion in most regions of interest and one must 
solve the evolution equations numerically starting from input distributions 
measured, for example, in deep inelastic scattering experiments [|HS|] . 

Suppose we wish to find the distribution function of W + bosons with 
helicity A, carrying a momentum fraction x, inside an electron neutrino, 
fw + lu e i x i Q 2 )- The lowest order expansion for this distribution function is 

fw+/u.frQ?)=Pw+/u.faQ a ) ( 5 ) 

where P w +, Ue (x, Q 2 ) is the splitting function which embodies the details of 
an elementary W + v e e~ vertex (see Appendix A). Keeping only this first 
term of the expansion is known as the effective vector boson approximation 
(EVBA) |25] which is the weak interaction analog of the effective photon ap- 



proximation. Weak splitting functions have been derived in many different 
contexts under various assumptions p6|, p7|, |28|, |29|. Different authors agree 



for terms proportional to ln(Q 2 /Mjy), the so-called "leading-log" contribu- 
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tions, but usually not further. Additional terms in the splitting functions can 
depend on the context of the derivation which can lead to process dependence 
and gauge dependence. If the EVBA is used in the context of, for example, 
Higgs production where Q 2 3> M^, the "leading-log" term in the splitting 
function is, in fact, the most important. However, our calculations require 
the use of more refined splitting functions, listed in Appendix A, which are 
slight variations of those of Johnson et a/.p8[. 



3.2 Naive Onset of Multiple Gauge Boson Emission 

The effective number of small-x partons associated with a high energy 
proton is closely related to the strength of the QCD coupling constant and the 
inability to solve the QCD evolution equations in a iterative manner starting 
from first principles. Likewise, the number of weak partons associated with 
a high energy fermion can suggest the suitability of keeping only the lowest 
order approximations to the weak parton distribution functions. Turning 
this argument around, the failure of the lowest order approximation to the 
weak parton distribution functions suggests an energy scale above which 
multiparticle processes become important. Neglecting, for the moment, the 
separate issue of restrictions which govern the feasibility of factorizing weak 
cross sections, it is amusing to consider two qualitative arguments leading to 
an estimate of such an energy scale. 

Consider a neutrino entering a reaction characterized by a center of mass 
energy ^/s. If, for example, the hard subprocess in the reaction involves 
the emission of a Z boson from the neutrino, the integral of the Z boson 
distribution function over all momentum fractions roughly corresponds to 
the number of Z bosons, N z / U (s), "inside" the neutrino 

N z /u(s)= I dxP z/u (x,s). (6) 

JMl/s 



The lower limit of integration accounts for the Z boson mass and it is assumed 
that ^fs is large enough to accommodate real Z boson emission. The precise 
kinematic limits are not of critical importance for our argument. In Fig. 1 
we plot Nz/ U { s ) which is summed over all boson polarizations although the 
transverse polarizations dominate the sum. It is relatively straightforward 
to show (see Appendix A) that as s — > oo, 

"'^toS^M*®' (7) 

which becomes of order unity at \fs ~ 10 8 GeV. We interpret this simple 
estimate as a scale above which the assumption of single Z boson emission 
is inappropriate; higher order corrections representing the possibility of mul- 
tiple emissions should then be considered. 

Had we used the "leading-log" approximation to the Z boson distribution 
function, we would have found 



47T COS 2 9„ 



Nz/u(s)* A _2,« ln 2 (,/Mi) (8) 



which is of order unity at y/s ~ 10 6 GeV. This lower scale reflects the well 
known shortcoming of the "leading-log" approximation for describing small- 
x fractions within the EVBAP3. For most applications discussed previously 



in the literature this small-x region is irrelevant but for us it is not. 

A second technique to estimate the onset of multiple gauge boson emission 
involves a consideration of the iterative solution to the evolution equations 
governing the parton distribution functions. In QCD this approach forces 
one to iterate to all orders because of the large coupling constant and the 
singular nature of splitting functions near the endpoints x = and x — 1. 
In the weak parton model the gauge boson masses avoid splitting function 
singularities and the small coupling suggests the validity of a series expansion. 
At asymptotically high energies the weak parton model should behave more 
like QCD in the sense that the gauge boson masses can be neglected. A crude 
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measure of when this occurs is to ask when the second order term in the weak 
parton model expansion of the distribution function becomes comparable to 
the first order term. 

The second order contribution to the Z distribution function in terms of 
first order splitting functions, depicted in Fig. 2, may be written as 

(2) r l - M z/ s dy 

fzL( x > s ) = — Pu/u(y,s) P z/u (x/y,ys) (9) 

Jx y 

where P u / U (y, s) = Pz/ui)- — V-, s ) corresponds to the emission of a Z boson 
by a neutrino. It is convenient to consider the ratio of the second order 
contribution to the first order (splitting function) contribution 

* J) = ^S' (10) 

The quantity R(x, s) is insensitive to x for small x, where the distribution 
functions are largest. The plot of R(x = lOMf /s, s) in Fig. 3 suggests that 
the second order correction becomes comparable to the first order term at 
v/i « 10 8 GeV. 

It is pleasing, though perhaps not surprising, to find agreement between 
the above two techniques which suggest the onset of multiple gauge boson 
emission. Nevertheless, these arguments must not be taken too literally 
since we have ignored the subtleties of whether weak factorization sanctions 
such calculations. Indeed, in our estimation of cr^j in the next section, we 
will encounter a context in which we must be considerate of factorization 
concerns. 

3.3 Parton Model Contributions to <7 t Q t ^ 

The elastic and inelastic Born cross sections for VeVn scattering at center 
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of mass energy y/s are given by 

l + 2z- 2z(l + z) log 

l + 2w-2w(l + u;)loi 



aelas [S) 4cos 4 ^Mf 



VeVji, , 



7ia 2 
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f w 
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where to = M^/s and z = M|/s. For a^ 1 ~ 30, M z = 91.17 GeV and M w = 
80 GeV these expressions for ct^'(s) and cr^ M (s) respectively asymptote to 
~ 70 pb and ~ 210 pb for y/s > 1 TeV. 

In order to estimate the energy at which o\^^ may begin to receive sig- 
nificant contributions from multiparticle production, consider the gauge in- 
variant set of diagrams of Fig. 4 corresponding to the exclusive production of 
an additional fermion pair in the reaction Vgp^ — > e~fi + ff. We shall focus on 
this process as being representative of a class of reactions at this order of per- 
turbation theory, though many other final states are certainly possible. Let 
us consider how these diagrams fit into a distribution function framework. 

In analogy with QCD, the electroweak factorization ansatz of Eq. [| asso- 
ciates the diagrams of Figs. 4a,b with a W + W~ fusion subprocess. Included 
in Fig. 4c are final state bremsstrahlung processes corresponding to fragmen- 
tation functions of the final leptons and initial state bremsstrahlung processes 
related to redefinitions of the neutrino structure functions. 

We can accommodate the 2 — ► 2 Born processes and the processes of 
Figs. 4 in the parton framework by summing only over W + W~ and z/ e z/ M 
subprocesses in the factorized expression for c^J in Eq. |3]. In addition, we 
only keep terms up to 0(a w ) in f w +,(x, Q 2 ) and f w -, p (x, Q 2 ) as in Eq. [|, 
and in 

f» e /v e &Q 2 ) = <Kl-x)+P ZK (l-x,Q 2 ), (13) 

/p M /p>,Q 2 ) = 6(l-x)+P z/Pli (l-x,Q 2 ), (14) 

where the delta functions reflect the initial conditions of the neutrinos. Con- 
sequently, in the factorization approximation, we use the above expansions 
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to obtain 

C>) = *?iZ n (s) + o*?( 8 ) + a^iZ-(s) + a brcm (s) + • • ■ (15) 

where the explicitly listed terms correspond to particular subprocesses which 
we will discuss in turn. 

The cross sections cr^°™(s) and cr^\ n (s) are given by Eqs. |Hi[l2| . The 
W + W~ fusion cross section is given by 

<&&-{*)= (16) 

J2 dx x dx 2 P w + /Ue (xi,s)P w - /i? (x2,x 1 s)a w + w - f f(x 1 x 2 s) 

x<d(x 1 x 2 s -4Jl4) 

where the sum extends over the polarization states of the on-shell W bosons. 
The subprocess cross sections for W^Wy — > // are straightforward to calcu- 
late and may be found in the literature [^TJ. The kinematic constraints of the 
parton model strictly limit the applicability of this particular contribution to 
fermion pair masses mjf > 2M^. The bremsstrahlung contributions cr brcm 
correspond to the emission of a near-mass-shell photon or Z which subse- 
quently decays into a // pair; since these only contribute to mjf < 2M\y 
(which falls outside the constraints of Eq. 16) we neglect them in this ap- 
proximation. 

In Fig. 5b we show cr^+°^_ summed over all lepton and quark species 
(except for t quarks) such that m^j > 2My/- The In s growth of this con- 
tribution to (Jtotai originates from the logarithms in the splitting functions 
corresponding to finding transverse gauge bosons inside neutrinos, in com- 
plete analogy with the QED expression of Eq. |2|. In the next section we 
compare this result with an exact calculation. 
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4 Multiparticle Production Beyond the Par- 
ton Model 

Fig. 5c gives the result of an exact calculation for fermion pair produc- 
tion for mjf > 2M\y which we performed using helicity techniques ||32|| . The 



helicity calculation incorporates interference effects between the diagrams of 
Fig. 4 and grows as Ins. As discussed in more detail in Refs. [^, Q the 



fact that the parton model behaviour disagrees with the exact calculation 
only by a power of In s and not by some power of s (which would be char- 
acteristic of a violation of gauge invariance in the EVBA) is somewhat of a 
success for the EVBA which neglects the effects of bremsstrahlung diagrams. 
Furthermore, closer inspection reveals that in the weak parton model the 
dominant contribution to // production originates from the threshold re- 
gion (rnw+w- — m ff > 2Mw) whereas the validity of the EVBA to W + W~ 
fusion requires that m w + w - 3> 2My/- It is at this level of agreement (or 
disagreement) that the weak parton model is useful for calculating low order 
contributions to cr^j. In any case, the growing contribution to cr^j from 
fermion pair masses above 2Mjy is not significant when compared with the 
much larger Born terms. 

Fortunately, we are not limited to using the weak parton model ex- 
clusively. We next consider particle production processes below the 2Mw 
threshold of the weak parton model. By integrating our exact 2 — » 4 cal- 
culation over all fermion pair masses we obtain a contribution to ct^j of 
0{a w /My / \a{s/My / )) 1 as shown in Fig. 5a, which is dominated by the Z 
boson pole. The same behavior is more easily deduced from the 2^3 pro- 
duction of an on-shell Z boson (see Fig. 6) which agrees with the 2^4 
calculation (except that the latter has an additional 0(o4) component cor- 
responding to the region \mjf — Mz\ > Tz)- 

In Fig. 7 we plot additional contributions to Gr^tai from the production 
of on-shell W ± and Higgs bosons (treating them as stable particles). At 

13 



■s/s ~ 10 7 GeV the sum of the 2^3 contributions is approximately half that 
°f ^da™ + °"in°i n — 290 pb. To be consistent, we should more appropriately 
compare the 2 — ^ 3 contributions with the radiatively corrected 2 — > 2 final 
states. To do this we consider the limit of the standard electroweak model 
without QED (i.e., sin6 w = 0) and take M w = M z = 91.17 GeV. In Fig. 8 
we show the sum of the analogous 2^3 contributions as well as an estimate 
of the of the radiatively corrected 2^2 contributions deduced from the 



work of Frankfurt and Sherman 33 



For the 2 — ^ 2 contribution to <J t ^ in Fig. 8 we take 

15vr</ A 

r 2 z W 



a 2 ^ 2 ~ a^(s) + aZns) - B(a - *o)^|m (£) (17) 



where <x^ M (s) and cr^ M (s) are given by Eqs. |TT| , |T2"| and the coefficient of the 
logarithm may be extracted, with the help of the optical theorem, from the 
leading log calculations of Ref. |3"3fl. As an additional approximation which 



fixes the nonleading logarithmic contribution to the radiative correction, we 
include the theta function and adjust the argument of the logarithm so that 
the radiative corrections only contribute above a center of mass energy t/sq. 
In this approximation the 2 — > 3 contributions surpass the 2-^2 contribu- 
tions to (Jtotai m the neighbourhood of 10 7 GeV for y/so in the range from 
M z to 10 TeV. 

The appearance of large logarithms in both the lowest order 2^3 pro- 
cesses and the next to leading order 2^2 processes raises concerns over the 
contributions of higher orders of perturbation theory. As a step towards es- 
timating the effects of higher orders we adopt a conjecture of Fadin, Kuraev 
and Lipatov of Ref. 0] whereby we modify our 2 — ^ 2 and 2^3 amplitudes 
in an attempt to sum leading logarithms to all orders. For the 2 — > 3 ampli- 
tudes, this corresponds to recognizing the dominance of the multiperipheral 
diagrams (e.g., Fig. 6a) and replacing the the vector meson propagators with 
those of reggeons. We should emphasize that unlike Ref. f|, it is not clear 
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that such replacements are sensible in our context. With this caveat in mind, 
we will show how our results are modified before discussing their significance. 
Consider the contribution to the 2 — > 3 amplitude from the diagram 
in Fig. 6a. Due to vector boson propagators this diagram is proportional 
to (qf — M 2 )~ l (q 2 — M 2 ) -1 with momentum transfers q\ = p Ue — p e - and 
<?2 = £>p m — Pft+- In the ansatz of Fadin, Kuraev and Lipatov these propagator 
factors are replaced by 

1 1 ( Sl /M 2 ) a(q2l) {S2/M 2 ) a{ql) 



q\ -M 2 q\- M 2 q\ - M 2 q\- M 2 



(18) 



where 



and 



Sl = (Pe~ +Pz) , S 2 = (p M + +p Z ) , (19) 

Q{q2) = J^fJ dkl M - M 2 )\(q - k)\ - M 2 ) ■ (20) 

Likewise, the 2^2 amplitudes are modified by an overall factor of (s/M 2 ) a( - q ) 
Incorporating these modifications in our helicity calculations leads to the 
cross sections shown in Fig. 9. The 2—^3 curves in Fig. 9 are extrapo- 
lated beyond 10 7 GeV since our unoptimized numerical integration becomes 
unreliable beyond that energy. Nevertheless, in this approximation we see 
that the 2^2 and 2^3 contributions to o\^ t ^ become comparable in the 
vicinity of 10 10 GeV. Again, these results are appropriate for sm.9 w = and 
M = M w = M z = 91.17 GeV. 

As mentioned above, the conjecture of Fadin, Kuraev and Lipatov de- 
serves scrutiny in our context. The implied modification of the 2 — > 2 ampli- 
tude most plausibly accounts only for the exchange of weak isospin T = 1 in 
the t — channel without manifestly accounting for the more important T = 
(pomeron) exchange. On the other hand, our estimate of the radiative cor- 
rections to 2 — ► 2 in Fig. 8 contains both T = and T = 1 contributions (but 
is only of first order in a w In s) . Hence it is perhaps safer to focus on the qual- 
itative nature of the higher order corrections suggested in Fig. 9. Namely, 
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that while higher order corrections tend to suppress the contributions of in- 
dividual channels, the energy scale at which 2 — > 3 reactions surpass 2 — > 2 
reactions is not extensively affected (to logarithmic accuracy). 

Though we have only explicitly dealt with 2 — > 2, 2 — > 3, and, to some 
extent, 2 — > 4 contributions to c^J, it is tempting to speculate beyond 
few body final states. Fig. 9 suggests that radiative corrections to 2 — > 3 
processes subdue the growth of those channels and may eventually lead to 
overall decreasing contributions. Meanwhile Pindor et al. || suggest a slowly 
growing cr^j which implies that 2 — > n (n > 3) channels must be opening 
up to compensate for the decreasing contributions from 2 — > 2 and 2 — > 3 
channels. This picture of multiparticle production is analogous to what is 
observed experimentally in QCD, albeit at a much different energy scale. 

In summary, we have shown that by considering the lowest order weak 
parton model and 2 — > 2 and 2 — > 3 contributions to <J t ^ that multiparticle 
production may become an important feature of weak interactions at energy 
scales of 10 7 — 10 9 GeV. In each case the appearance of large logarithms 
plays an important role yet we suspect that the relevant energy scale is not 
an artifact of a finite order of perturbation theory. 
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Appendix A 

For completeness, we include in this appendix explicit expressions for 
the vector boson splitting functions used in the text. The following formulae 



are essentially those of Johnson et al. from Ref. p8[ with some minor but 
notable modifications. 

Consider a fermion entering a reaction with a body A at center of mass 
energy y/s (see Fig. Al). The splitting function formalism addresses processes 
which correspond to the emission of a vector boson by the fermion so the 
vector boson is considered as an initial state particle in a subprocess with 
the body A. If the lagrangian coupling between vector bosons and fermions 



is vl/r^W, where 



9r1h 



7s 



fit 7m' 



75 



2 u "* 2 ' 

then the longitudinal vector boson splitting function is 

"2(1 -x)(r + x) 



P^ 



V L /f 



X,S 



'9l + 9r 



16ir 2 J x{x — r) 

2r(2 + r-x) 



(21) 



(22) 



(x — r)'' 



In 



(1 + r — x)x^ 



where r = My/s accounts for the vector boson mass My. The spin factor 
F s = 2 for neutrinos (and F s = 1 otherwise) compensates for the fact that 
the formulae of Ref. [^8| are spin averaged assuming the initial fermion has 
two helicity states. 

The splitting functions correspond to the emission of a positive or negative 
helicity boson are, respectively, 



v+/f 



x,s) 



9 2 L hi{x,r) +g 2 R h 2 (x,r) 



(23) 



Pv-/f{x,s) = F s g 2 R h 1 (x,r)+g 2 L h 2 (x,r) 



(24) 
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where 
hx(x, r) 



167T 2 



^1 - x)(2 + r-x) 



[x — r) 

[1 + r — x) 



— (x — r) ln(x) 



(25) 



x + r \ f(l + r — x)x 



(x — r) 5 



h 2 (x,r) 



16vr s 



(1 — x)(2 + r — x) x + r /(1+r — x)x^ 
(1 + r — x)(x — r) (x — r) 2 y r 



(26) 



Eqs. p2~ , p5 ,26 are slightly different from the corresponding expressions of 
Ref. ||28|| . Our use of an on-shell flux factor has the effect of replacing an 



overall factor of x in the equations of Ref. |28j with a overall factor of x — r. 
A consequence of this modification is that the distribution functions vanish 
at x = r. 

The quantity Ny/f{s) very roughly corresponds to the probability of find- 
ing a vector boson (of any polarization) "inside" a fermion. It is defined as 
the integral over x of the vector boson distribution functions, which to lowest 
order, is given by 



iV, 



v//l 



dx 



p v L /f(x, s) + Pv+/f{x, s) + Pv-/f{x, s) 



Using the above expressions for the splitting functions results in 



iV, 



v/f\ 



' 9l + 9 2 R 

167T 2 



,2 7, r 2 r 19 
In r H — In r + — — or -\ 



(27) 



(28) 
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Figure Captions 



Figure 1. Integral of the lowest order Z boson distribution function for a 
neutrino. 

Figure 2. Kinematics relevant to second order contribution to the Z boson 
distribution function for a neutrino, s is the invariant mass squared 
of the original vA system. One Z boson carries away a momentum 
fraction 1 — y whereas the Z boson entering the hard subprocess carries 
a fraction x of the original neutrino momentum. 

Figure 3. Ratio of second order contribution to lowest order Z boson dis- 
tribution function for a neutrino. 

Figure 4. Gauge invariant set of diagrams for the production of an addi- 
tional fermion pair (where / ^ //, e). 

Figure 5. (a): a{v e v^ — > e~fi + Z) calculated exactly at tree level; (b) aivep^ - 
e~[x + ff) in effective vector boson approximation (ruff > 2M\y)', (c) 
oivep^ —* e~n + ff) calculated exactly at tree level (mff > 2Mw)', 

Figure 6. Diagrams for the production of an on-shell Z boson, (a) Dom- 
inant multiperipheral contribution (b) Subdominant bremsstrahlung 
diagrams. 

Figure 7. 2 — > 3 contributions to C7 t „ t ^ in the standard model with Mz = 
91.17 GeV, M w = 80 GeV and M H = 500 GeV. Cross sections for final 
states containing a Higgs boson (a(y e p ll H') + cr(e~fx + H)), W boson 
(a(u e fi + W~) + a{e^u^W + )), Z boson (cr(e~{i + Z) + a{u e u^Z)) and 
their sum. 

Figure 8. Contributions to o\^^ in the limit sin# w = 0, My/ = Mz = 
Mn — 91.17 GeV. Sum of 2 — >• 2 contributions contains estimated 
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lowest order radiative corrections which are assumed to vanish below 
below ^/s = 1 TeV. 2-^3 curve is sum of exact tree level calculation 
for states containing final state Z, W or H bosons. 

Figure 9. Contributions to u^tai m the limit sin 9 W = 0, M w = M z = 
Mn — 91.17 GeV with amplitudes modified by the Fadin, Kuraev and 
Lipatov (FKL) conjecture described in the text, (a) Sum of 2 — > 2 
Born contributions, (b) Sum of FKL modified 2 — > 2 contributions, (c) 
Sum of 2 — ► 3 Born contributions, (d) Sum of FKL modified 2^3 
contributions. The dashed portions of (c) and (d) are extrapolations. 

Figure Al. Diagram corresponding the process fA — > f'X in which the 
initial fermion / emits a gauge boson V which, in the effective vector 
boson approximation, is considered to be on-shell in the subprocess 
AV^X. 
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